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ABSTRACT 

We present deep photometry in the B,V and I filters from CTIO/MOSAIC for about 270.000 stars in the Fornax dwarf Spheroidal 
galaxy, out to a radius of r f //as0.8 degrees. By combining the accurately calibrated photometry with the spectroscopic metallicity 
distributions of individual Red Giant Branch stars we obtain the detailed star formation and chemical evolution history of Fornax. 
Fornax is dominated by intermediate age (1-10 Gyr) stellar populations, but also includes ancient (10-14 Gyr), and young (<1 Gyr) 
stars. We show that Fornax displays a radial age gradient, with younger, more metal-rich populations dominating the central region. 
This confirms results from previous works. Within an elliptical radius of 0.8 degrees, or 1.9 kpc from the centre, a total mass in stars 
of 4.3xl0 7 M Q was formed, from the earliest times until 250 Myr ago. 

Using the detailed star formation history, age estimates are determined for individual stars on the upper RGB, for which spectroscopic 
abundances are available, giving an age-metallicity relation of the Fornax dSph from individual stars. This shows that the average 
metallicity of Fornax went up rapidly from [Fe/H]<-2.5 dex to [Fe/H]=-1.5 dex between 8-12 Gyr ago, after which a more gradual 
enrichment resulted in a narrow, well-defined sequence which reaches [Fe/H]a;-0.8 dex, «3 Gyr ago. 

These ages also allow us to measure the build-up of chemical elements as a function of time, and thus determine detailed timescales 
for the evolution of individual chemical elements. A rapid decrease in [Mg/Fe] is seen for the stars with [Fe/H]>-1.5 dex, with a clear 
trend in age. 

Key words. Galaxies: dwarf - Galaxies: evolution - Galaxies: stellar content - Galaxies: Local Group - Stars: C-M diagrams 



1. Introduction 

Due to their close proximity to the Milky Way (MW), dwarf 
spheroidal galaxies (dSph) can easily be resolved into individual 
stars, and studied in detail. Once thought to be the simplest 
systems in the Local Group (LG), the dSphs hav e been found 
to dis play complicated evolutionary histories (e.g.. lTolstov et al.l 
° l2009h . 

Fornax is a well studied galaxy, discovered by Shapley in 1938, 
and is one of the most luminous and massiv e dSphs in the Local 
Group, second only to the Sagittarius dwarf (Ma teoll998l) . It has 
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an absolute magnitude of My=— 13.3, and a total (dynamical) 
mass of « 1.6x10 s M , from modelling o f kinematic data out 
to a radius of 1 .9 k pc dWalker et al.l 12006: Klimentowski et aD 
2007; iLokasI 12009). Fornax is located at a distance of 138+8 
kpc ((m-M)v=20.84±0.04), determined using RR Lyrae stars, 
in good agreement with other measurements using the in- 
frared tip of the Red Giant B ranch (RGB) method and the 
Horizontal Branch (HB) lev el dGreco et alJ l2005t iRizzi et al.l 
120071 : IPietrzvriski et alJl2~0()9). The estimated red dening toward 
Forna x, from extinction maps, is E(B-V)=0.03 ( Schlege l et al.l 
fl998T) . 

Fornax can easily be resolved into individual stars, but given 
the substantial size on the sky (r t ;d a i=71±4 arcmin or 2.85+0.16 
kpc), wide-field imaging is need ed to study the spatial distri- 
bution of the resolve d stars (e.g.. llrwin & H atzidimitriou 1995; 
Battaglia et al. 2006). The first s tructural study of Fornax was 
made with photographic plates by Hodge ( 1961b), revealing that 
the ellipticity of Fornax isophotes increases with distance from 
the centre. Furthermore, an asymmetr y was found with a peak 
densi t y offset from the centra l position (Ide Vaucouleurs &" Abies 
11968b iHodge & Smithlfl974l) . Unusually for LG dwarf galaxies , 
several globular clusters (GC) were also found (Shapley 1938; 
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lHodgefl961al) . around the Fornax dSph 
CCD studies were later conducted in the Forna x dSph 
down to ever decreasing brightness levels ([Eskridgel 1 1988 



alpha, iron-peak and heavy-elements. 



Demers et al.lll994Hlrwin & H atzidimitriou 1995; Walcher et 
20031) . The presence of old stellar populations (>10 Gyr) 
was confirmed by the dis covery of several RR Lyrae 
stars ( Bersier & Wood 2002), and the presence of an ex- 
tended HB~JSOTlan^^Ll|2003). Photometric studies have also 
revealed a large number of carbon stars, ind icating a substantial 
intermediate age (1-10 Gyr) popul ation dAaronson & Mouldl 
1980, 11985b lAzzopardi etalj 1 19991) . In the central regions 
of Fornax, young Main Sequence stars (<1 Gyr) were first 
identified by Buon anno et al.l (1 19851) . which led to the discovery 
of a stella r population gradie nt of the young and intermediate 
age stars (Ste tson et al.lll998l) . Photometric surveys of Fornax 
have found stellar over-densities and "shells", some of which 
are beli eved to be the result of a recent encounter with a sm aller 
system dColeman et al.ll200 4, 2005b; lOlszewski et al.|[2006l) . 
The Star Formation History (SFH) of Fornax has 

been determined using de ep C olour-Magnitude 

dStetson et all Il998t iBuonanno et all 



Diagrams (CMDs ) 



1999: ISaviane et alJl2000b IGallart et alJl2005h . Using wide-field 
CMDs Hstetson et alj ( 19981) determined the age of the young 
M S population to be 10 0-200 Myr, with similar results found 
bv lSaviane et alJ (|2000). Deep photometry in small fields close 
to the centre of Fornax indicate that star formation started « 12 
Gyr ago, continuing almost to the pres ent day, with a young 
episode of star for mation 1-2 Gyr ago dBuonanno et alJ ll999: 
Gallart et al. 2005). Furthermore, hints were found of separate 
Sub-Giant Branches, indicati ng that star formation may have 

occurred in discrete episodes dBuonanno et al.ll 19991). 

A dee p photometric survey was made by [Coleman & de Jongl 
d2008h . covering an area of 5.25 deg 2 centred on Fornax, down 
to a 50% completeness limit of B=23.0 and R=23.5. The SFH 
was determined at different radii from the centre, showing 
that Fornax experienced a complex evolutionary history with 
numerous epochs of star formation. A significant population 
gradient was found with radius, with ancient stars being present 
at all radii, but more recent star formation concentrated to the 
centre. The star formation is found to have gradually decreased 
from the earliest times, until a sudden episode of strong star 
formation occurred s4 Gyr ago. 

Numerous spectroscopic studies have also been performed in 
Fornax, of individual stars on the upper RGB. Medium resolu- 
tion Ca ii triplet spectroscopy has provided [ Fe/H] measurement s 
for increas ing numbers of s tars [33 star s dTolstov et all [2001). 



2001 



117 stars dPont et al 
675 stars dKirbv et al 
tailed Metallicity Distribution 
which shows that the dominant 



870 stars dBattaglia e taij 12008b . 
2010)]. These studies provide a de- 
Function (MDF) of Fornax, 
population is relatively 



metal-rich ([Fe/H]x:-0.9 dex), but that stars as metal-po or as 
[Fe/H]<— 3.7 dex are also present dTafelmever et alJ l2010). By 
analyzing the spatial distribution of stars in different age ranges, 
it became clear that the ancient stars display a more extended 
spatial distribution (with colder kinematics) than the intermedi- 
ate ag e stars, and the latter of the younger stars (Battaglia et al. 

High resolution (HR) spectroscopic studies of individual RGB 
stars in Fornax have been ca rried out within the central 0.3 de- 
grees of Fornax, for 3 stars (Tol stoy et al J 12003 : Shetrone et al] 
l2003l) . 9 stars in three globular clusters (iLetarte et all f2006) 
and more recently 81 stars dLetarte et al.l 1201 Oh using 
VLT/FLAMES. These studies have revealed the complex 
abundance pattern for individual stars in Fornax, including 




Fig. 1: Coverage of the photometric and spectroscopic observations 
across the Fornax dwarf spheroidal galaxy. The big squares denote the 
coverage of the CTIO 4m/MOSAIC fields, while the smaller squares 
show the CTIO 0.9m fields. The open (blue) circ les show the stars ob- 
served in the VLT/FLA MES Can triplet survey dBattaglia et alJl2008t 
Starkenbur g et all 120101) . for which only [Fe/H] is deter mined. The 
solid (red) dots mark the stars observed us ing medium dKirby et all 



120101) and high resolution spectroscopy dLetarte et al] 120101) . for 
which both [Fe/H] as well as a-element abundances are determined. 
T he dashed (green) ell ipse is the tidal radius of Fornax, as given 
byHattaelia et aQ i2006h . 



3V e 

IT 



The observed properties of the Fornax dSph have been 
modelled several times, in simulations using different 
techniques, making Fornax o ne of the best studied dSphs 
in the Local Group ( e.g.. iLanfranchi & Matteuccl [2003 



Mas hchenko et al l 120081: iRevaz et al.l 120091; iKirbv et alJ 12011 
iRevaz & Jablonka 2012]). The chemo-dyn amical Smoothed- 
Particl e Hydrodynamics (SPH) code used bv lRevaz & Jablonkal 
d2012l) correctly matches the observed metal-rich MDF of 
Fornax, and shows a narrow [Mg/Fe] distribution, consistent 
with trends derived from HR spectroscopy. Cosmological 
simulations have reproduced the stellar parameters of a model 
galaxy in good agreem ent with the Fornax dSp h, including a 
globular cluster system (Mas hchenko et aH l2008). Furthermore, 
chemical evolution models of the Fornax dSph have been able 
to reprod uce the observed MDF and alpha-element dis tribu- 
tion (e.g JLanfranchi & Matteuccill2003l: IKirbv et al.ll201 ll) . 
In this work, we present new accurately calibrated, wide-field 
photometry of the Fornax dSph, going down to the Main 
Sequence Turn-Off (MSTO), carried out with MOSAIC on the 
CTIO 4m/Blanco telescope. Using the accurate photometry, the 
SFH is determined over a larg e area of the Fornax dSph, using 
CMD synthesis me thods ( e.g.. iTosi et alJJ 199 it i Tolstov & Saha 
1996 ; IGallart et alj|1996ah lDolphinlll997t lAparicio et al.l ll997: 
Ide Boer et al] 2012 ). Spectroscopic metallicities (from Can 
triplet spectroscopy) are directly used in combination with the 
photometry, to provide additional c onstraints on t h e age of the 
stellar populations, as desribed in Ide Boer et alj (|2012). The 
detailed SFH is used to determine the probability distribution 
function for age of individual stars, giving age estimates for 
stars on the upper RGB. These ages are linked to observed 
spectroscopic metallicities and abundances, giving an accurate 
Age-Metallicity Relation (AMR) and the timescale of chemical 
enrichment in the Fornax dSph. 
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The paper is structured as follows: in section [2] we present our 
photometric and spectroscopic observations. Section [3]presents 
the CMDs of Fornax and an analysis of the spatial variations 
of stellar populations within Fornax. In section [4] we describe 
the specifics of the method used to obtain the SFH of Fornax. 
The analysis of the detailed, spatially resolved SFH is given in 
section [5] and the chemical evolution timescales derived from 
individual RGB stars in section [6] Finally, section [7] discusses 
the results obtained from the SFH and chemical evolution 
timescales. 



2. Data 

2.1. Photometry 

Deep optical photometry of the Fornax dSph in the B, V and I fil- 
ters was obtained using the CTIO 4-m MOSAIC II camera over 
9 nights in October 2008 and November 2009. Deep photom- 
etry for the outer fields was obtained as part of observing pro- 
posal 2008B-0397 (PI M. Mateo) supplemented by deep point- 
ings covering the centre of Fornax as part of observing proposal 
2009B-0157 (PI A. Saha). Our observing strategy was to obtain 
several non-dithered, long (600s) exposures for each pointing, 
which were stacked together to obtain the deepest photometry 
possible. Short (10s, 90s) exposures were also obtained, in order 
to sample the bright stars that are saturated in the deep images. 
In order to ensure accurate photometric calibration of the data 
set, several fields were observed in Fornax using the 0.9m CTIO 
telescope, under photometric conditions. Further more, observa- 
tions were also made of Landolt standard fields dLandom [l992. 
[2007b covering a range of different airmass and colour. 
The positions of the observed fields from the 4m telescope are 
shown in Figure Q] The spatial coverage of the B and V filters is 
complete for r cW <0.8 degrees, while the I filters is complete for 
r f /;<0.4 degrees. The position, exposure time, airmass and see- 
ing conditions of the Fornax fields are given in Appendix lAl 
The reduction and accurat e calibration of this d ataset follows the 
steps described in detail in de Boe r et all (1201 ll) . An accurate as- 
trometric solution was determined for each pointing, after which 
the different exposures of each field were aligned and stacked 
together to obtain a single, deep imag e. Photometry was ca rried 
out on all images using DoPHOT dSchechter et al.lll993h . An 
accurate photometric calibration solution was determined using 
the standard star observations, depending on airmass, colour and 
brightness. All fields were calibrated and subsequently placed 
on the same photometric scale using the overlapping regions 
between different fields. Finally, the different fields were com- 
bined to obtain a single, carefully calibrated photometric cata- 
log, which is given in table 2. 

2.2. Spectroscopy 

Medium and high resolution spectroscopic observations are 
available for individual stars on the RGB in Fornax, giving mea- 
surements of [Fe/H] and the detailed abundance patterns of indi- 
vidual stars. 

Medium resolution (R~6500) Ca n triplet spectroscopy is avail- 
able for 870 individual RGB stars in the Fornax dSph, from 
VLT/FLAMES observations dBattaglia et al.ll2006l 120081) wit h 
the Ca ii triplet [Fe/H] calibration from Starkenbu rg et al.l yOlO). 
These observations provide [Fe/H] measurements for stars out to 
a radius of 1.3 degrees from the centre of the Fornax dSph, and 
include a range in metallicities from -3.0<[Fe/H] <-0.5 dex. 
Furthermore, medium resolution spectroscopy by Kirb v et all 



(12010 ) provides [Fe/H] and a-element abundances for 675 stars 
within the central 0.3 degrees of Fornax, with a metallicity range 
of -2.8<[Fe/H]<-0.1 dex. From this dataset, only those stars 
with an uncertainty in [Fe/H] and [Mg/Fe] lower than 0.2 dex 
are used, to select only the brightest stars with a high signal-to- 
noise. This brings their sample to a total of 72 stars to be consid- 
ered in the analysis. 
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Fig. 2: Comparison between [ Fe/H] and [Mg/Fe] from HR and medium 
resolu tion spectroscopy, from Letarte et al. ( 2010) (HR), Battagli a et all 
(2008) (CaT) and high S/N stars of iKirbv et all d2010h (Kirbv) us- 
ing overlapping stars, a) Comparison of the metallicity of HR spec- 
troscopy with high S/N Kirby (red) and CaT (blue) stars, as a function 
of metallicity. b) The difference between medium resolution spectro- 
scopic datasets of Kirby and CaT. c) Comparison between the [Mg/Fe] 
determined using HR and medium resolution spectroscopy. 



Furthermore, HR spectroscopy (from VLT/FLAMES) is also 
available forr e //<0. 3 degrees, for 81 individual RGB stars in the 
Fornax dSph dLetarte et al.l 120 10), providing [Fe/H] as well as 
c-elements (O, Mg, Ca, Si, Ti) and r- and s-process elements (Y, 
La, Ba, Eu, Nd). For [a/Fe] obtained from these observations, we 
assume [a/Fe] =( [Mg/Fe] +[Ca/Fe] +[Ti/Fe])/3. The HR spec- 
troscopy covers a range in metallicity from -2.6<[Fe/H]<-0.6 
dex, for stars on the upper RGB. 

The spatial coverage of the spectroscopic observations is shown 
in FigureQ] for both the DART HR and Ca n triplet spectroscopy 
as well as the Kirby spectroscopy. 

To determine if the DART and Kirby spectroscopic catalogs 
are consistent, we compare [Fe/H] and [Mg/Fe] abundances for 
overlapping stars. Figure [2] s hows tha t the metallicities derived 
for Ca ii triplet datasets Batta glia et al.l d2008l);ISta rkenbur g et al.l 
(2010) and medium resolution spectroscopy IKirbv et al.l (2010) 
are consistent over the metallicity range probed by the overlap- 
ping stars. Furthermore, Figure [3J; shows th at the [Mg/Fe] de- 
rived for the highest S/N measurements from lKirbv et all (1201 01) 
ar e consistent with tho se determined from the HR spectroscopy 
of lLetarte et all (1201 0l) . 

3. The Colour-Magnitude Diagrams 

The full, calibrated photometry catalog discussed in Section |2~T1 
contains photometry in the Fornax dSph in the B,V and I filters 
for radii out to r e ;;«1.5 degrees. Due to the spatial coverage of 
the catalog, the B and V filters are fully sampled for radii r f «<0.8 
degrees, while the I filter is complete only for r e //<0.4 degrees. 
Figure|3]shows a (V, B-V) CMD of the inner r c //<0.4 degrees of 
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Fig. 3: Fornax (V, B-V) CMD of the central r f/ /<0.4 degrees with boxes identifying the different stellar populations. BHB: Blue Horizontal 
Branch, RHB: Red Horizontal Branch, RGB: Red Giant Branch, MSTO: Main Sequence Turn-Off, BP: Blue Plume, BL: Blue Loop, RC: Red 
Clump, AGBb: AGB bump. A box providing an estimate of the Milky Way foreground contamination is also shown (FG). Error bars represent the 
photometric errors at each magnitude level. 



the Fornax dSph, marked with the different features that can be 
identified in the CMD. 

The CMD of Fornax displays a wide variety of evolutionary 
features, consistent with an extended period of star formation. 
The presence of blue RGB stars and an old MSTO indicates 
ancient stars (>10 Gyr), which is consistent with the presence 
of Blue Horizontal Branch (BHB) stars. A well populated Red 
Horizontal Branch (RHB) and Red Clump (RC) indicate that the 
main population of stars in Fornax has intermediate age (1-10 
Gyr), further supported by the presence of a strong red RGB 
population. Clear signs of recent star formation (<2 Gyr) are 
also found in Figure [3] such as the presence of Blue Loop (BL) 
stars and a strong Blue Plume (BP) population. 
To determine the photometric completeness of the CMD, a large 
number of artificial star test simula tions have been carri ed out, 
similar to the method described in de Boe r et al. I (|20T1 . A set 
of artificial stars was generated with parameters within the range 
0.25<age<15 Gyr, -2.5<[Fe/H]<-0.30 dex, -0.2<[a/Fe]<0.60 
dex. Stars were distributed randomly across the six MOSAIC 
pointings in Fornax, in 3 different filters. To avoid changing the 
crowding properties in the images no more than 5% of the total 
observed stars were ever injected as artificial stars at one time. 
A total of 400 images in each of the two inner fields, and 200 
images in each of the four outer fields were created, containing 
5000 artificial stars each. This resulted in 4000 images, contain- 



ing a total of 8 million artificial stars spread across the full area 
of the Fornax dSph. 
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Fig. 4: The 50% completeness level in the B (solid blue line), V (green 
dashed line) and I (red dotted line) filters for increasing elliptical ra- 
dius (r e n) from the centre of the Fornax dSph. The vertical (black) line 
indicates the radius up to which observations in the I filter are com- 
pletely sampled. 
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0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 2 

B-V B-V B-V 

Fig. 5: (V, B-V) Colour-Magnitude Diagrams of the Fornax dwarf spheroidal, shown for increasing elliptical radius. Error bars showing the 
average photometric error and a 50% completeness line (red) are also included. The solid blue line shows a reference isochrone corresponding to 
the RGB of the dominant population ([Fe/H]=-1.00 dex, [a/Fe]=0.00 dex, Age=4 Gyr). 




0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 

VI VI VI 

Fig. 6: (I, V-I) Colour-Magnitude Diagrams of the Fornax dSph for increasing elliptical radius. The error bars show the average photometric error 
at each magnitude level, and the red line indicates the 50% completeness level. The solid blue line shows a reference isochrone corresponding to 
the dominant RGB population in Fornax. 

Figure |4] shows the variation of the 50% completeness level with The completeness levels for the B and V filters are comparable 
distance from the centre of Fornax, in the three different filters, at each radius, while the 50% completeness level in the I filter 
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is systematically brighter due to the increased sky brightness in 
I. However, the higher intrinsic brightness of stars in the I band 
compensates for the brighter completeness level. The complete- 
ness level in the I filter is shown only for r e //<0.4 degrees, due to 
the lack of I band observations in the outer parts of Fornax. The 
centre of Fornax is less complete at a fixed magnitude limit than 
the outer regions, due to the large amount of crowding in the 
central region, which means fewer stars are unambiguously de- 
tected, and placing the 50% completeness at brighter magnitude 
levels. 



3.1. Spatial variations 

Using the full catalog we can construct CMDs at different radii 
from the centre of the Fornax dSph, which can be used to study 
the radial behaviour of different evolutionary features. Figures[5] 
and|6]show carefully calibrated (V, B-V) and (I, V-I) CMDs for 
the Fornax dSph at different elliptical radii (i e u). Error bars show 
the photometric errors at each magnitude level, and the (red) line 
indicates the 50% completeness level at each elliptical radius. 
A reference isochrone tracing the dominant RGB of Fornax is 
also shown, in blue. No de-reddening has been applied to the 
CMDs in Figures [5] and [6] Instead, models and isochrones used 
in the CMD analysis are reddened using the extinction coeffi- 
cient adopted for Fornax (Schlegel et al ll998h . 
Figures [5] and [6] show that the stellar populations in Fornax 
change with radius from the centre. Ancient stars (>10 Gyr) are 
visible at all radii, indicated by the presence of a blue RGB and 
old MSTO stars. Additionally, the outskirts of Fornax show a 
well populated BHB population, which is heavily contaminated 
in the inner parts by the presence of young stars. 
The red RGB and RC populations are dominant in the centre of 
Fornax, but disappear from the CMDs with increasing elliptical 
radius. This is a clear sign that the intermediate age stars are 
more centrally concentrated than the old stars. 
The bright BP stars are confined to the innermost radial bin in 
Figures [5] and [6] The outskirts of Fornax do show a popula- 
tion of blue stars extending from the old MSTO, but these stars 
display a different slope than the BP stars observed in the cen- 
tral region. This feature could be due to the presence of young 
stars, but is likely caused by th e presence of Blue Straggler 
Stars (BSS) dMapelli et alJ l2009h . The bright BP populations 
disappear rapidly at larger distances from the centre, which in- 
dicates a radially changing age and/or metallicity of stars in 
Fornax. 

This is also seen in the radial distribution of the BP, RGB and 
MSTO populations, as given in Figure [7] The young (BP) pop- 
ulation is most centrally concentrated, while the stars on the 
RGB (which traces both the intermediate and old populations) 
are more extended. Finally, the old MSTO stars show an ex- 
tended radial distribution, with stars present at all radii. To study 
the stellar populations of Fornax, and their trends, in more detail 
we analyse the different CMD features and search for variations 
with elliptical radius. 



3.1.1. The Red Giant Branch 

The RGB of Fornax displays a wide spread in colour, as seen 
in Figures [5] and [6] This is in dicative of a spread in metallic- 
ity and/or age, as first found bv lStetson et al.l (Il998l) . To investi- 
gate the properties of RGB stars in detail we use the medium 
resol ution Can triplet spectr oscopy of individual RGB stars 
from Starkenburg et al. (2010) (see Section l2~2l i. Figure[8]shows 



the RGB of Fornax overlaid with the [Fe/H] metallicities of in- 
dividual RGB stars where Can spectroscopy is available. Due to 
the higher crowding in the central region a brighter magnitude 
limit was adopted for the spectroscopic observations in the in- 
nermost radial bin in Figure [8] 




Fig. 7: Normalised radial distribution of the BP, RGB and MSTO pop- 
ulations (shown in Figure[3} of the Fornax dSph. 



The inner regions of Fornax show stars with metallicities rang- 
ing from -2.8<[Fe/H] <-0.2 dex. The strong red RGB compo- 
nent seen in Figures |5]and|6]shows a metallicity of [Fe/H]«-1.0 
dex, and is well traced by an isochrone with an age of 4 Gyr. 
For increasing radius, Figure [8] shows that the most metal-rich 
populations ([Fe/H]>-0.7 dex) rapidly disappear. The dominant 
RGB population also diminishes with radius from the centre, 
albeit at a slower pace. Conversely, the metal-poor, old popu- 
lations ([Fe/H]<-1.7 dex) remain present at all radii. This ob- 
served effect with distance from the centre indicates the presence 
of a metallicity gradient with radius in Fornax. 

3.1.2. The Main Sequence Turn-Offs 

The MSTO region of the Fornax dSph (see Figure [5} shows the 
presence of stellar populations with a wide range of ages, from 
very young (BP) to very old (red MSTO). Previous studies of this 
region have found clear evi dence of a population gradient with 
radiu s in the Fornax dSp h dStetson et alJll998t FBuonann o et al.l 
1999: lGallartet al.ll2005h . Here, we present deep homogeneous 
photometry of the MSTO region covering a large fraction of the 
galaxy out to r e //=0.8 degrees. 

Figure [9] shows the MSTO region of the Fornax dSph 
at different radii from the centre. Reference isochrones 
have been overlaid on the CMDs, representative of the 



present in Fornax ( 


Battadia et al J20081 IStarkenburs et al J20 1 Ot 


Letarte et al. 2010 


). The ages of the isochrones have been se- 



lected to fit the colours of RGB stars of the corresponding metal- 
licity or the observed distribution of young MS stars. Figure [9] 
clearly shows that the youngest populations are confined to the 
centre of Fornax, while the older populations are more extended. 
The youngest BP populations disappear rapidly as we go out- 
ward from the Fornax centre. 

For increasing elliptical radius the dominant intermediate pop- 
ulation also diminishes, as becomes clear when comparing the 
area underneath the (red) metal-rich reference isochrone at dif- 
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Fig. 8: The RGB region of the Fornax dSph at different elliptical rad ii, x M . The (black) dots s how the photometry, while the larger (coloured) 
circles indicate the [Fe/H] measurements of individual RGB stars from Starkenbu re et al.l d2010h . 



ferent distance from the centre. Conversely, the oldest popula- 
tions traced by the (blue) metal-poor isochrone remain present 
at each elliptical radius. A random selection of equal numbers 
of stars in each panel shows that this effect is not due to the de- 
creased total number of stars at different elliptical radii. These 
trends are consistent with a radial gradient in both metallicity 
and age in the Fornax dSph, spanning a wide range in age. 

4. SFH Method 

To determine the SFH of the Fornax dSph, we compare the 
observed CMDs with a grid of synthetic CMDs using Hess 
diagrams (plots of the density of observed stars) of stars 



... — ™ — 

thetic CMD method (e.g., Tosi et al. 


1991; Tolstov&Saha 


1 19961 iGallart et al.lll996al iDolphin 1997 


: lAparicio et al. 


1997). 



Uniquely, Talos simultaneously takes into account observations 
in all available photometric filters as well as the spectroscopic 
MDF , to obtain a well constrained SFH, as discussed in detail 
in lde Boer eTaIl(l2012l) . 

First, a set of synthetic ideal CMDs and MD Fs is generated usin g 
the Dartmouth Stellar Evolution Database (iDotter et al. 2008), 
for population bins covering the range in age, metallicity and a- 
element abundances. To simulate observational conditions in the 
synthetic CMDs used to determine the Fornax SFH, we make 
use of the artificial star test simulations described in [3] This 
approach is the only way to take into account the complex ef- 
fects that go into the simulation of observational biases, such as 



colour-dependence of the completeness level and the asymme- 
try of the photomet ric errors of stars at faint magnitudes (e.g. 
IGallart et alJ U996b). The results from the artificial star tests are 
used to include observational errors in synthetic CMDs in a sta- 
tistical manner, so they can be compared directly to the observed 
CMDs, to obtain the best matching SFH. 

Using the synthetic CMDs, model MDFs are generated by ex- 
tracting only the stars corresponding to the same magnitude 
range as the spectroscopic sample. The stars are binned in metal- 
licity and convolved with a Gaussian profile to simulate obser- 
vational errors. 

The difference between model and observed CMD and MDF 
is mi nimized according to a Poisson equivalent of x 1 {Dolphin 
120021) . The uncertainties on the SFH are obtained by determin- 
ing the SFH using a se ries of different CMD and p arameter grid- 
dings, as described by Apari cio & Hida lgo (2009j). The average 
of all solutions is adopted as the final SFH, with the standard de- 
viation as error bars. In this way, the best matching star forma- 
tion rate (SFR) is obtained for each population bin (with age and 
metallicity) considered in the SFH. These values are then pro- 
jected onto one axis to obtain the SFR as a function of age (Star 
Formation History) or metallicity (Chemical Evolution History). 
Talos has been extensively tested, and has been shown to accu- 
rately regrochac£fte_r£operties of both real and synthetic test 
data ( de Boer et al.ll2012l) . Before applying Talos to our photo- 
metric and spectroscopic data sets of the Fornax dSph, we need 
to take into account the specifics of the SFH determination of 
this galaxy. 
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Fig. 9: The distribution of stars in the MSTO region of the Fornax dSph, for different elliptical radii. Three isochrones are overlaid, representative 
of the metal-rich ([Fe/H]=-1.00 dex, [o'/Fe]=0.00 dex, age=4 Gyr, red) and metal-poor populations ([Fe/H]=-2.45 dex, [ff/Fe]=0.40 dex, age=14 
Gyr, blue) on the RGB, as well as the young population of the Fornax dSph ([Fe/H]=-0.30 dex, [or/Fe]=-0.20 dex, age=0.3 Gyr, green dotted). 



4.1. Parameter space 

The limits of populations adopted in the SFH are determined 
by considering all available information from the current data 
sets (see Section[2]i as well as the literature. 
For the limits on metallicity, we consider the spectroscopic 
MDF from Can triplet and HR spectroscopy. On the metal- 
poor end, the [Fe/H] is limited to [Fe/H]=-2.5 dex, by the 
availabilit y of metallicity in the Dartmouth Stellar Evolution 
Database dDotter et al]l20 08). However, this is not a problem, 
since the MDF shows that almost no st ars in the Fornax spectro - 
scopic sample have [Fe/H]<-2.5 dex (IStarkenburg et al.l l2010). 
On the metal-rich end, the spectroscopic MDF shows that no 
stars with [Fe/H]>-0.5 dex are present on the RGB. However, 
higher metallicities may be present in young Main Sequence 
stars, which have no corresponding RGB. From previous deriva- 
tions of the SFH, no stars have been found m ore metal-rich than 
rFe/Hl=-0. 1 dex (Cole man & de Jongl2 008). which we adopt as 
the metallicity limit in our SFH. A binsize of 0.2 dex is assumed 
for [Fe/H], which is similar to the average observed uncertainty 
on [Fe/H]. 

In order to take into account the c-element abundances, we 
adopt a linear [Fe/H] vs [q/Fe] relation, b ased on the HR spec- 
troscopic abundances (Letart eet al.ll20l6T) . A range in [q/Fe] is 
assumed for each bin in [Fe/H], which reproduces the trends in 
the HR data. 

The choice of possible ages to adopt in the SFH solu- 



tion is motivated by previous derivations of the SFH of the 
Fornax dSph. Previous work has shown that Fornax con- 
tains stars covering a ll ages, from ancient (>10 Gyr) to very 
young (<1 Gyr) age s ( Buona nno etalll 1 9991: iGallart et al JI2005I: 
IColeman & de Jong||2008l) . This is also consistent with the ob- 
served CMDs, which show a prominent blue plume, caused by 
stars with ages less than 1 Gyr. The lower limit in age is set 
by the adopted isochrone set, which does not contain isochrones 
with ages<0.25 Gyr. Assuming a maximum age of 14 Gyr, for 
the age of the Universe, we consider therefore a range of ages 
between 0.25 and 14 Gyr old. A binsize of 1 Gyr is used for ages 
>5 Gyr, 0.5 Gyr for ages between 1-5 Gyr and 0.25 Gyr for 
ages <1 Gyr, to take into account th e different age sensitivity of 
young MSTOs (Hidalgo^eL^|2011|). 



4.2. Age resolution 

The age resolution of the SFH solution determined by Talos is 
important to understand the limitations of our derived SFH of 
Fornax. The ability to resolve episodes of star formation at dif- 
ferent ages depends mainly on the photometric depth of the ob- 
served photometry. Furthermore, the method of determining the 
uncertainties of the solution (see Section |4| results in a smooth- 
ing of the SFH, which gives limits to the resolution of the SFH at 
different ages. The age resolution of the solutions is determined 
by the ability to recover the S FH of a set of synth etic populations 
at different input ages (e.g. JHidalgo etall2011h . 
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Fig. 10: Input and recovered SFHs of a series of synthetic short (10 
Myr) episodes of star formation at young (a), intermediate (b) and 
old (c) ages. The black, solid histogram shows the input SFH, given 
the adopted binning of the solution. The red histograms show the re- 
covered SFH, as well as the fit of a Gaussian distribution as the dashed 
line. The mean (pi) and variance (<x) of the fitted Gaussian distribution 
are also listed. 



A set of seven synthetic short episodes of star formation (with 
a duration of 10 Myr) was generated, with different input ages 
covering the age range adopted for Fornax (0.6, 2.2, 4.2, 6.5, 
8.5,10.5 and 12.5 Gyr). The metallicity distribution of the syn- 
thetic episodes is chosen to match the observed MDF of Fornax. 
Using artificial star test results, observational conditions were 
simulated corresponding to the central region of Fornax. Given 
that the crowding is highest in this part of the galaxy (and the re- 



sulting photometry the least deep), this will result in upper limits 
to the age resolution. The SFH for the synthetic episodes is re- 
covered by fitting the photometry in three filters (B,V and I), 
simultaneously with a synthetic 50% complete MDF with simi- 
lar photometric depth as the Ca n triplet spectroscopy. 
The recovered and input SFH for the synthetic episodes is shown 
in Figure[10l for young (a), intermediate (b) and old (c) star for- 
mation episodes. The recovered SFH is fit by a Gaussian distri- 
bution, from which we determine the age of the central peak (//) 
as well as the variance cr, which determines the resolution of the 
recovered episode. Figure [10] shows that the central peak is re- 
covered at the correct age for young ages (<7 Gyr), given the 
binning adopted for the solutions. For older ages the peak of the 
synthetic episode is recovered at slightly too old ages. 
The Gaussian fits show that the age resolution of the recovered 
episodes is significantly better at young ages than at old ages. 
For the young ages with small age resolution, the star formation 
is confined mostly to the central bin, while for the old episodes 
the star formation is spread out over multiple bins. The episodes 
are recovered with a resolution of «0.1 Gyr at an age of 0.6 
Gyr, ^0.4 Gyr at an age of 2.2 Gyr, ^0.8 Gyr at an age of 4.2 
Gyr ^1.5 Gyr at an age of 6.5 Gyr, »1.8 Gyr at an age of 8.5 
Gyr, w2.3 Gyr at an age of 10.5 Gyr and ^2.5 Gyr at an age of 
12.5 Gyr, which is consistent with values between 17 and 22% 
of the adopted age. 

5. The Star Formation History of the Fornax dSph 

Having described our method for determining the SFH 
of the Fornax dSph in Section |4] we now apply Talos 
to our photo metric and spectros c opic data sets of the 
Fornax dSph dBattariia et al] l2008t IStarkenburg et al l 120101: 
iLetarte et al]|2010tlKirbv et al.ll2010h . 

The SFH is derived using the photometric CMDs in combina- 
tion with the spectroscopic MDF, to constrain the final SFH as 
described in Ide Boer et al.1 d2012l) . Within the central 0.4 de- 
grees, photometric information in three filters (B,V and I) is 
combined with the available HR and Can triplet spectroscopy. 
For radii r e ;;>0.4 degrees, the SFH is determined using just the 
(V, B-V) CMD combined with Can triplet spectroscopy. The 
spectroscopic MDF is used to constrain the SFH solution for 
-2.5<[Fe/H]<-0.6 dex. CMD analysis shows that more metal- 
rich stars are too young to appear on the RGB. Therefore, these 
populations are only constrained by the photometry. An inves- 
tigation of the MDF at different photometric depth has shown 
that no luminosity bias is present in the spectroscopic sample of 
Fornax stars. 

Figure QT| presents the final SFH and Chemical Evolution 
History (CEH) of the Fornax dSph out to a radius of r e //=0.8 de- 
grees, with error bars indicating the uncertainty on the SFR as a 
result of different CMD and parameter griddings (as described in 
Section©. The SFH and CEH display the rate of star formation 
at different ages and metallicities over the range of each bin, in 
units of solar mass per year or dex respectively. The total mass 
in stars formed in each bin can be determined by multiplying the 
star formation rates by the range in age or metallicity of the bin. 
The SFH shows that star formation is clearly present at all 
ages, from as old as 14 Gyr to as young as 0.25 Gyr. Most 
of the star formation takes place at intermediate ages, between 
1-10 Gyr, consistent with earlier determinations of the Fornax 
SFH (G allart et al.ll2005t IColeman & de Jong|l2008l) . The CEH 
in Figure ITTb shows that the dominant population displays 
metallicities between -1.5<[Fe/H]<-0.7 dex, but that signifi- 
cant levels of star formation are also present at [Fe/H]x:-2.4 dex 
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Fig. 11: The overall a) Star Formation History and b) Chemical Evolution History (CEH) of the Fornax dSph, out to an elliptical radius of 
r c i/=0.8 degrees. The SFH and CEH have been determined using all available photometric information, incombination with observed spectroscopic 
metallicities. 



and [Fe/H]^-0.2dex. 

Figure [12] shows the overall SFH of Fornax, divided into 
young (a), intermediate (b) and old (c) age ranges, with the cor- 
responding CEH. The youngest star formation in Fornax (<1 
Gyr) displays declining SFR for younger ages, with a metal-rich 
CEH showing a peak at [Fe/H]»-0.2 dex. The ages of these pop- 
ulations are determined with a high accuracy (see Section I4.2l i. 
indicating that the change in SFH corresponds to a real drop in 
SFR. 

The intermediate age star formation shows a double peaked SFH 
with a declining SFR for ages <6 Gyr. The resolution of the 
SFH at these ages (<2 Gyr) is good enough to distinguish be- 
tween both peaks, indicating the presence of distinct episodes 
of star formation. Furthermore, the age extent of the younger 
episode (<6 Gyr) is inconsistent with a single short burst of 
star formation and indicates a smooth change in star forma- 
tion activity over several Gyrs. The metallicities of the popu- 
lations with ages between 2-5 Gyr (solid) display a range be- 
tween -1.5<[Fe/H]<-0.5 dex with a peak at [Fe/H]»-1.0 dex. 
The slightly older intermediate age star formation (5-10 Gyr, 
dashed) is dominated by more metal-poor stars with a distri- 
bution between -2.1<[Fe/H]<-0.7 dex, showing two peaks at 
[Fe/H]«-1.8and-1.2dex. 

Finally, the oldest ages (10-14 Gyr) display a declining SFR 
with increasing age, and a metal-poor CE H with a similar to - 
tal range as observed in the Sculptor dSph dde Boer et al .1 l20l2h . 
However, the SFR at the oldest ages is not constrained well due 
to the depth of the observed CMDs, resulting in a smoothing 
of the underlying actual SFH and removing any observable fea- 
tures. 

To check the distribution of stars in the CMD we compare the 
best matching observed and synthetic CMDs of the Fornax dSph 
in Figure [13] The synthe tic CMD has been generated using the 
Dartmouth isochrone set dDotter et al.l l2008). and therefore does 
not include a RC or HB. The overall distribution of stars in the 
observed CMD is well matched in the synthetic CMD. Stars on 
the BP and MSTO occupy the same range in colour and bright- 
ness as in the observed CMD. Furthermore, the colour range and 
slope of the RGB is a good match to the observations, as shown 
by comparison to the reference lines. 

The total number of stars in the synthetic CMD is consistent with 
the observations to within a few percent. Additionally, the rela- 
tive fraction of stars occupying the RGB and BP phase is similar 



to that of the observed stars. This shows that the total mass in 
stars in the Fornax dSph is well matched, and distributed across 
the CMD with the right numbers. Therefore, we feel confident 
that the SFH constitutes an accurate representation of the overall 
behaviour of the stellar populations in the Fornax dSph. 
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Fig. 12: The overall SFH of the Fornax dSph, divided into young (a), 
intermediate (b) and old (c) age ranges. The corresponding CEH is also 
shown for each age range. For intermediate age stars, younger (solid) 
and older (dashed) stars are separated to show the effect of age on the 
CEH. 



Using the best matching SFH we can show the age and metallic - 
ity corresponding to different CMD features. Figure [T4l shows 
the synthetic (V, B-V) CMD of Fornax generated using the 
Dartmouth isochrones, colour coded with age and metallicity for 
each individual star. The effect of different stellar populations on 
the MSTO region of the CMD is clearly seen from Figure [14] 
with younger Main Sequence stars occupying bluer and brighter 
positions on the BP. The CMD shows that the majority of stars 
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is of intermediate age, but that lower levels of old star formation 
are also present. The position of different stellar populations on 
the RGB shows that the oldest stars occupy the blue side of the 
RGB, while progressively younger stars are found toward the red 
side of the RGB. However, small numbers of young (<4 Gyr) 
stars are also found on the blue side of the RGB, after evolving 
off the young MS. 

The effect of metallicity on the positions of stellar populations 
is shown in Figure [T4b . The RGB is sensitive to [Fe/H] rather 
than age, and shows a clear metallicity pattern across the RGB, 
which provides a good match to the observed metallicity pattern 
seen in Figure[8] Stars on the BP are mostly metal-rich, although 
metal-poor stars are also present. This could be due to the lack of 
spectroscopy of BP stars, which means the metallicity of these 
populations is not constrained from the MDF fitting. 

5. 1. The effect of different isochrone sets 

Figure Q3] shows the SFH solution of the central region 
of Fornax (r f //<0.1 94 degrees) obtained using three differ- 
ent isochrone sets dDotter et all l2Q08h iPietrinferni et all 12004 
lYi et al.l l2()01). The central region was chosen since the range 
of populations is greatest here, allowing a comparison of deter- 
mined ages over a large timespan. Figure [T5l shows that the so- 
lution is very similar in all cases, and consistent within the error 
bars. The overall behavior of the SFH is recovered in all solu- 
tions, with a peak at *4 Gyr, slowly declining SFR toward older 
ages, and a flat distribution for ages >10 Gyr. The similarity 
between different isochrone sets is likely due to the dominance 
of the young Main Sequence and MSTO in the solution (where 
isochrone agree reasonably well), while the contribution of the 
RGB (which is where isochrone sets disagree most) is taken into 
account solely through the MDF. Therefore, the shape of the 
RGB is not fit in the solution, but rather the number of stars 
occupying the upper RGB (as determined by the isochrone lu- 
minosity function). The good agreement of the SFH solutions 
gives confidence to our final SFH of Fornax. 

5.2. Spatial variations in the SFH 

The coverage of the photometric and spectroscopic information 
extends to a radius of r e //=0.8 degrees. Therefore, it is possible to 
derive the SFH at different radii from the centre, which allows us 
to determine the presence of variations in the SFH with radius, 
as seen in Section [3~Tl 

The total radial extent of the photometry and spectroscopy is di- 
vided into 5 annuli, each containing a similar number of stars ob- 
served in the B and V filters, since these filters are fully complete 
out to greater radii from the centre. For each of the annuli, the 
SFH and CEH are determined by fitting all available photomet- 
ric and spectroscopic information, in the same way as discussed 
in Section [5] For the inner three annuli, photometric informa- 
tion in all three available filters (B, V and I) is used, together 
with the spectroscopic MDF. The SFH in the outer two annuli 
is determined using just the (V, B-V) CMD, together with the 
spectroscopic MDF from Ca n triplet spectroscopy. 
Figure[l6]shows the best matching SFH and CEH in each of the 
five annuli of the Fornax dSph. Figure [To*h shows that a clear 
gradient with age is present in Fornax, with the peak of star for- 
mation shifting to progressively older ages for increasing radius. 
The CEH shown in Figure [T6b shows that the distribution of 
metallicities changes from peaked in the inner parts to a broad 
distribution in the outermost bin. 



The young star formation (age<l Gyr) is present only within the 
centre of Fornax, consistent with the observed CMDs. This is 
also seen in the CEH, where significant star formation at high 
metallicities ([Fe/H]>-0.5 dex) is found only in the central part 
of Fornax. Very low levels of young star formation are seen in 
the outer parts, which are likely due to fitting the probable BSS 
population in the outer parts as a young stellar population. 
The intermediate age population of Fornax displays a clear gra- 
dient in age with radius from the centre. The peak of star forma- 
tion shifts from an age of ^8-10 Gyr in the outermost annulus 
to a peak at age^4 Gyr in the central bin. The CEH of Fornax 
changes from a broad distribution in the outer regions to a peak 
at [Fe/H] «- 1.0 dex in the centre, in good agreement with the 
dominant RGB population seen in Figures [5] and [6] Given the 
resolution of the SFH solution at these ages (see Section 14.21 1. 
this radial change in the SFH corresponds to a smoothly varying 
SFR with radius in Fornax, inconsistent with distinct short bursts 
of star formation. 

Old star formation (age>10 Gyr) is present at all radii in Fornax, 
with a flat distribution as a result of the poor SFH resolution. 
The SFR starts to increase gradually in all re gions at an age 
of &10 Gyr, consistent w ith previous studies (Buonanno et alJ 
[19991: lGallart"eTaT]l2005h . 

The observed and synthetic MDFs inferred from the SFH are 
shown in Figure [T7J for all five radial annuli. The synthetic 
MDFs are consistent with the observed spectroscopic MDFs 
within the errors, as expected, given that the MDF is used as an 
input in the SFH determination. In the outer parts of Fornax, the 
metal-poor component is slightly over-estimated in the synthetic 
MDF, while the inner parts show that the metal-rich populations 
are slightly under-estimated. The MDF on the upper RGB shows 
the same trend with radius as the CEH in Figure [T6b . chang- 
ing from a peaked distribution in the innermost bin to a more 
broad distribution in annulus 5, consistent with the diminishing 
strength of the dominant intermediate age component. 
The total mass in stars formed in our SFH of the Fornax dSph 
is 4.3xl0 7 M within 0.8 degree or 1.9 kpc. The core radius 
of the mass distribution resulting from the SFH solution is 
r c =0.29±0.03 degrees, determined using a Sersic profile fit. This 
is consistent within the error s with the core radiu s derived from 
the observed density profile dBattaglia et alj |2006). 



6. The timescale for chemical evolution of the 
Fornax dSph 

Using the SFH presented in Figures QT| and [16] it is possible 
to determine the age probability function of individual stars 
on the RGB. For each observed star, all stars in the syn- 
thetic CMD (generated using the SFH) with the same mag- 
nitude (in all filters) and metallicity within the observed un- 
certainties are selected. The ages of these synthetic stars build 
up a probability distribution function for age, which is con- 
sidered to be re presentative o f the a ge of the observed star, 
as described in Ide Boer et alJ ([2012). Assuming a Gaussian 
distribution, the mean of the distribution function is adopted 
as the age of the observed star, with the standard deviation 
used as an error bar. Using this method, we determine accu- 
rate age estimat es for all available samples of spectroscopic 
stars in Fornax (Battaglia et alJ I2.Q08I : IStark enburg e t alJ 120101 : 
iLetarte et alfeolOUKirbv et al.ll2010h . 
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Fig. 13: The a) observed and b) synthetic Hess diagrams for the Fornax dSph, out to r ( , / ;=0.8 degrees. The colours represent the number of stars 
in each bin, on a logarithmic scale. To allow a better comparison between both Hess diagrams, (black) reference lines are also shown. The BHB, 
RHB and RC are not present in the synthetic Hess diagram, since they are not modelled in the isochrone set used. 



a) age 




1 



8 g 
CD 

0! 

6 5" 




Fig. 14: The synthetic CMD of the Fornax dSph inferred from the best matching SFH as determined using the Dartmouth isochrone set. The 
colours represent a) age and b) [Fe/H] for each individual synthetic star. 
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Fig. 15: Comparison between the a) SFH and b) CEH of the centre of Fornax (r e/ /<0. 194 degrees), determined using different isochrone sets. The 
SFH solution is shown for Darthmouth Stellar Evolution Database isochrones (red, solid), the BaSTI/Teramo isochrones (blue, dotted) and the 
Yonsei-Yale isochrones (green, dashed). 



6. 1 . Age-metallicity relation 

Using the accurate ages for stars with spectroscopic metal- 
licities, we determine the detailed AMR of the Fornax 
dSph (Figure [18]). The results are consistent with previous deter- 



minations of the AMR dPont et al]l20"ol: lattaglia et al. 2006), 
but determine the age of stars with greater accuracy, reveal- 
ing the detailed effects of different populations. The AMR 
shows that the metallicity in Fornax increased rapidly from 
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Fig. 16: The SFH and CEH for 5 annuli of Fomax containing roughl y equal number of stars, with the radial extent indicated in each panel. For 
comparison, the core radius of the Fornax dSph is 0.23 degrees (Irwin & Hatzidimitriou 1995). The SFH for annuli 1-3 is obtained by fitting all 
available photometry (in the B, V and I filters) and spectroscopy. For annuli 4 and 5 only the V,B-V CMD is used together with the spectroscopy 
to constrain the SFH. 



[Fe/H]<-2.5 dex to [Fe/H] =-1.5 dex between 8-12 Gyr ago. 
However, the poor SFH resolution results in a large scatter in 
age, removing any features at the oldest ages. 
At intermediate ages (3-8 Gyr), the interstellar medium is grad- 
ually enriched during several gigayears, resulting in a narrow, 
well-defined AMR which reaches [Fe/H]x:-0.8 dex at ^3 Gyr. 
Furthermore, from Figure [18] hints are seen that the slope 
of the AMR may become steeper at even younge r ages (<3 
Gyr), corresponding to a phase of rapid enrichment (Pon t et alJ 
2004). However, more spectroscopic observations, in particu- 
lar of young BP stars, are needed to unambiguously detect this 
change in slope. The young, metal-poor stars in the lower left of 
Figure [18] most likely correspond to AGB stars which are fit as 



if they were RGB stars (due to the absence of the AGB phase in 
the adopted isochrone set), leading to a too young age estimate 
for these stars. 

A comparison to the AMR of the Sculptor dSph shows that the 
slope in Fornax at the oldest ages is steeper than in Sculptor, 
which indicates that metal enrichment occurred with a different 
rate. After the initial phase of metal enrichment, the slope of the 
AMR is comparable in both galaxies, but Sculptor only shows 
stars with old (>7 Gyr) ages. 

Figure [19] shows the spatially resolved AMR of Fornax, in each 
of the five annuli described in Section l5T2l The shape of the AMR 
is consistent in all annuli, with a similar slope at intermediate 
ages. At the oldest ages, more scatter is seen in the inner two 
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annuli than in the outer three annuli. The presence of the age 
gradient is visible in Figure [19] as an increased number of stars 
with younger ages in the AMR for more centrally concentrated 
regions. 



covered by the metal-rich stars in Figure l20h indicates that the 
"knee" will likely be spread over a large metallicity range, which 
may be a result of the extended SFH and multiple populations of 
Fornax. 



6.2. Chemical abundances 



7. Discussion 



The accurate age estimates can also be coupled to the stars for 
which other elements are determined, to derive the evolution of 
individual elements as a function of time. In this way, we can 
directly measure the detailed timescale for chemical evolution 
in Fornax. Figure l20h shows the [Mg/Fe] measurements for in- 
dividual RGB stars as a function of [Fe/H], with the age of the 
individual stars colour coded, and Figurel20b shows [Mg/Fe] di- 
rectly as a function of age. 




Fig. 17: Observed (red, dashed) and synthetic (green, solid) histogram 
MDFs from the SFH analysis in each of the five annuli in Fornax. 



Figure [20] shows the effect of different stellar populations on 
[Mg/Fe] in Fornax, and the rate with which this abundance 
changes at different times. Metal-poor stars are not well sampled 
in Figure [20] but the most metal-poor stars are consistent with 
the abundances observed for the MW halo (Lemasle et al., in 
prep). Conversely, the metal-rich stars are different from the MW 
halo, which indicates that Fornax experienced a different evolu- 
tionary history than the systems that may have for med the MW 
halo, as seen in other Local Group dwarf galaxies (Tolsto y et al.l 
2009). The metal-rich stars in Figure l20h display rapidly de- 
creasing [Mg/Fe] abundances, with a clear trend in [Fe/H] and 
age. Hints are seen of a changing slope in the [Mg/Fe] distribu- 
tion for populations of different ages. 

Connecting the metal-poor and metal-rich populations observed 
in the [Mg/Fe] distribution in Figure [18] a "knee" may be ex- 
pected to occur at [Fe/H] «— 1.5 dex, at an age of 7-10 Gyr. 
The "knee" in [Mg/Fe] marks the time at whic h SNe la start 
to contribute to the [Fe/H] content of a galaxy dTinslev 19791; 
Matteucci & Brocatol[l990l: iGilmore & Wvsd 1 199 it iMatteuccil 



We have presented the detailed SFH and CEH of the Fornax 
dSph out to r e //=0.8 degrees, obtained using a combination of 
deep, multi-colour photometry and spectroscopic metallicities. 
The obtained, spatially resolved SFH (see Figures [TTI and [Tol l 
shows features similar to previous SFHs (e.g. lGallart etaiT2 005: 
IColeman & de Jongl2008l) . However, the SFH presented here re- 
solves the stellar ages with greater accuracy, due to the inclusion 
of the spectroscopic MDF In particular, the age of the dominant 
episode of star formation is determined more precisely, consis- 
tent with the position of the RGB in the observed CMD (see 
Section l3XTT i. 




2003). However, due to the limited sampling of intermediate 
metallicity stars in Fornax, the "knee" is not visible in the data 
presented in Figure [20] Additionally,the large range in [Fe/H] 



Age (Gyr) 

Fig. 18: The Age-Metallicity Relation of stars on the upper RGB in 
the Fornax dSph, incorporating the full SFH a nd MDF information . 
Medium and hig h resolution spectroscopy from iLetarte et alj (2010); 
iKirbv et al] d2010h is shown a s blue points, while Can triplet spec- 
troscopy from Batt aglia et al] d2008l) is shown in red. The AMR 
of the Sculptor dSph, fr om HR spectroscopy, is shown as green 
points (de Boer et all2012h . 



The SFH of Fornax shows little variation in SFR at old ages (>10 
Gyr), and lacks a dominant episode of star f ormation at the old- 
est age s, contrary to what was observed bv lColeman & de Jong! 
(2008). The absence of a peak in the SFR at the start of star 
formation is likely an effect of the age resolution of the SFH, 
smoothing out any features at the oldest ages. To determine if the 
SFR peaks at the oldest age, deeper photometric data is needed. 
Fornax displays an extended SFH, dominated by star forma- 
tion at intermediate age (1-10 Gyr). The reason for the oc- 
currence and strength of this intermediate age episode is not 
clear. One possibility is a merger with a gas-rich companion, 
which triggered a strong episode of star formation. However, this 
gas would need to be pre-enriched at a relatively high metallic- 
ity ([Fe/H]«-1.5 dex). Another scenario is that initially expelled 
gas remained bound to the system for an extended period of time, 
and fell back onto the galaxy, leading to an episod e of star for- 
matio n similar to what is predicted by simulations dRevaz et alj 
2009). The fall-back of expelled gas could be caused by an in- 
teraction with the MW (Piateketal. 2007), similar to models 
proposed for the Carina dSph dPasetto et al.ll201 II) . Furthermore, 



14 



T.J.L. de Boer et al.: Star Formation History of the Fornax dSph 



the re-accretion of gas can also be a result of the orbit shape of 
Fornax, compressing and re-gathering expelled gas at apogalac- 
ticon, allowing for rep eated bursts of star formation over the 
course of several orbits (Nic hols et alJl2012h . 



gas may be responsible for the substructures, and may have trig- 
gered centrally concentrated star formation with a different spa- 
tial distribution. To investigate this link, the detailed kinematics 
and ages of young stars in Fornax need to be determined. 
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Fig. 19: The Age-Metallicity Relation of stars on the upper RGB in 
Fornax, in each of the five annuli, from medium an d high resolu- 
tion spectroscopy dLetarte et all2 010: Kirhv et al l 201 Oh as blue points, 
and Can triplet spectroscopy (Battaglia et al. 2008) as red points. 



Figure Q~2] shows that star formation was still taking place in the 
very centre of Fornax a few 100 Myr ago, resulting in a visible 
young population (see Figure [5]), which ranges in age between 
~2 Gyr and 0.25 Gyr. The age of the youngest stars may be 
younger than can be fit using the adopted isochrone set, which 
would lead to an incorrect estimate of the SFR in the youngest 
bin (0.25-0.50 Gyr). However, the level of star formation is 
expected to be low, not significantly influencing the results at 
young ages in Figure [T6l 

The young stars of Fornax display a different spatial distri- 
bution with respect to the old stars. Furthermore, the ages of 
these young stars ove rlap with those of the shell features found 
in the Fornax dSph (IColeman et alJ 12004. 2005a). Therefore, 
a link may be present b etween these features, as suggested 
bv IColeman et al.l d2005al) . The infall of a smaller system or of 




7.1. Comparison to the Sculptor dSph 

It is useful to compare the results deriv ed for Fornax wit h 
our previous study of the Sculptor dSph (Ide Boer et al ] l2012h . 
Sculptor is a very different galaxy, which experienced a much 
shorter evolution history, leading to a different SFH and abun- 
dance pattern. However, the same types of measurements were 
obtained in Fornax and Sculptor, allowing a detailed comparison 
between the properties of both galaxies. 

A comparison between the AMR of Fornax and Sculptor (see 
Figure [T8l ) suggests that the initial metal enrichment in Fornax 
occurred on a shorter timescale than in Sculptor, leading to 
a steeper slope at earlier times. A possible explanation for 
this could be the greater mass of Fornax, which allows it 
to retain more of its gas after SNe explosions, leading to 
more rapid e nrichment, similar to wh at is predicted by simu- 
lations (e.g., iRevaz &"ja blonka 2012). This rapid metal enrich- 
ment could also be responsible for the different position of the 
[Mg/Fe] "knee" in Fornax, co mpared to the observed "knee" in 
Sculptor (Ide Boer et al.ll2012l) . 

The dominant stellar component of the Fornax dSph is an inter- 
mediate age population. This is very different to Sculptor, which 
formed the bulk of its stars before «10 Gyr. This lack of inter- 
mediate age star formation may be due to its lower (dynamical) 
mass, which did not allow it to retain as much of its gas as Fornax 
and prevented it from continuing to form stars. 

8. Summary 

We have presented deep B, V and I CMDs of the Fornax dSph, 
out to an elliptical radius of 0.8 degrees, containing a total of 
270.000 stars. The CMDs show that Fornax displays a radial gra- 
dient with age and metallicity, with more metal-rich, younger 
stars forming more toward the centre. Fornax is dominated by 
intermediate age (1-10 Gyr) star formation, and shows a strong 
dominant RGB linked to a stellar population with an age of »4 
Gyr,at[Fe/H]«-1.0dex. 

Combining the deep CMDs with the spectroscopic MDF of indi- 
vidual RGB stars we have derived the detailed SFH of Fornax, at 
different radii from the centre. The spatially resolved SFH con- 
firms and quantifies the radial age and metallicity gradient, with 
a peak of star formation that shifts to older ages for increasing 
radius. From the SFH we determine that the total mass in stars 
formed in Fornax is 4.3xl0 7 M within an elliptical radius of 
0.8 degrees, or 1.9 kpc. 

Using the SFH, we determined age estimates for individual RGB 
stars, for which spectroscopic abundances are available. This al- 
lows us to study, for the first time, the evolution of different el- 
ements in Fornax, directly as a function of time. Using the in- 
dividual age estimates, we determined the detailed AMR of in- 
dividual stars and timescale for evolution of a-elements of the 
Fornax dSph. The AMR shows that Fornax experienced an ini- 
tial period of star formation, which enriched the metallicity from 
[Fe/H]<-2.5 dex to [Fe/H]=-1.5 dex from »12 to 8 Gyr ago. 
Subsequent gradual metal enrichment resulted in a narrow, well- 
defined AMR during the next 5 Gyr that reaches [Fe/H]«-0.8 
dex at *3 Gyr. Hints are seen in the AMR of a second change in 
slope at young ages (1-3 Gyr). 

The detailed timescale of chemical enrichment (see Figure l20b 
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[Fe/H] (dex) Age (Gyr) 

Fig. 2 0: a) [Mg/Fe] measurements for the medium and high resolution spectroscopic samples of RGB stars from iLetarte et alj|2010l : [Kirby et al. 
2010) (coloured filled circles). The colours represent the age in Gyr, derived from the SFH. Stars in the Milky Way are shown for comparison (small 
grey points), b) [Mg/Fe] is plotted directly against age for the same sample. 



shows that the abundances of old, metal-poor stars are con- 
sistent with the MW, while the abundances of more metal- 
rich stars show a departure from MW abundances. Metal- 
rich ([Fe/H]»-1.3 dex) stars show a rapid decrease of [Mg/Fe], 
with a clear trend with age. Connecting both features predicts 
the presence of a knee in the alpha-element distribution at 
[Fe/H]*- 1.5 dex, at 7-10 Gyr. Hints are seen of a different slope 
in the [Mg/Fe] distribution for different populations, which could 
be a result of the different strength of previous stellar popula- 
tions, leading to a different rate of chemical enrichment. 
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Appendix A: Observing log 



Table A.l: List of observed fields in the Fornax dSph, with the 4m CTIO Blanco telescope during two observing runs. Information is given about 
exposure time, airmass and seeing conditions as determined on image. 



Date 


Field 


RA 


DEC 


Filter 


exp time 


seeing 


airmass 






J2000 


J2000 








// 




2009 Nov 19 


Fnxcl 


02:39:51.84 


-34:47:48.96 


B 






1.2 


1.08-1.25 


2009 Nov 19 








V 


3600, 90, 


10 


1.2-1.4 


1.30-1.70 


2009 Nov 22 








I 


4800, 90, 


10 


1.0-1.5 


1.25-1.70 


2009 Nov 21 


Fnxc2 


02:39:51.84 


-34:13:48.96 


B 


3600, 90, 


10 


0.9-1.2 


1.08-1.26 


2009 Nov 21 








V 


3600, 90, 


10 


1.1-1.5 


1.30-1.70 


2009 Nov 23 








I 


4800, 90, 


10 


1.0-1.3 


1.20-1.68 


2008 Oct 4 


FnxNE 


02:42:07.28 


-34:03:34.90 


B 


6600, 90, 


10 


0.9-1.7 


1.07-1.55 


2008 Oct 4 








V 


9000, 90, 


10 


1.0-1.5 


1.09-1.74 


2008 Oct 6 


FnxNW 


02:37:25.24 


-34:11:44.40 


B 


4200, 90, 


10 


1.3-1.6 


1.00-1.21 


2008 Oct 6 








V 


7800, 90, 


10 


1.3-1.6 


1.00-1.15 


2008 Oct 5-6 


FnxSE 


02:42:28.77 


-34:42:00.10 


B 


5400, 90, 


10 


1.1-1.7 


1.00-1.70 


2008 Oct 5-6 








V 


7200, 90, 


10 


0.9-1.6 


1.00-1.80 


2008 Oct 4-5 


FnxSW 


02:37:52.57 


-34:50:12.30 


B 


4200, 90, 


10 


1.0-1.4 


1.08-1.43 


2008 Oct 4-5 








V 


7200, 90, 


10 


0.9-1.3 


1.06-1.60 


2008 Oct 5-6 


FnxIW 


02:39:51.84 


-34:47:48.96 


B 


1800 




1.2-1.4 


1.09-1.20 


2008 Oct 5-6 








V 


6000 




1.3-1.6 


1.06-1.30 
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